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ABSTRACT 

We use the very large Millennium Simulation of the concordance ACDM cosmogony to cali- 
brate the bias and error distribution of Timing Argument estimators of the masses of the Local 
Group and of the Milky Way. From a large number of isolated spiral-spiral pairs similar to the 
Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to 
true mass to be a factor of 1.8, while the 5% and 95% points of the distribution of this ratio are 
separated by a factor of 5.7. Here we define true mass as the sum of the "virial masses i\'jf200 
of the two dominant galaxies. For current best values of the distance and approach velocity of 
Andromeda this leads to a median likelihood estimate of the true mass of the Local Group of 
5.27 X 10^2 or log A/lg/M© = 12.72, with an interquartile range of [12.58, 12.83] and 
a 5% to 95% range of [12.26, 13.01]. Thus a 95% lower confidence limit on the true mass of 
the Local Group is 1.81 x 10^^ M0. A timing estimate of the Milky Way's mass based on the 
large recession velocity observed for the distant satellite Leo I works equally well, although 
with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of 
2.43 X 10^2 with a 95% lower confidence Hmit of 0.80 x lO^^ j^q. 

Key words: Galaxy: formation - galaxies: Local Group - galaxies: kinematics and dynamics 
- dark matter. 



1 INTRODUCTION 

During the 1970's it became generally accepted that most, per- 
haps all, galaxies are surrounded by extended distributions of 
dark matter, so-called dark ha los ( lEinasto. Kaasik & Sa ail 1 1974 
lOstriker. Peebles & Yahiilll974t) . These were soon understood to 
play an essential role in driv ing the formation and clustering of 
galaxies jWhite & R"ee3l 19781) . With the introduction of the Cold 
Dark Matter (CDM) paradigm, these ideas took more concrete 
form, allowing qua ntitative predictions to b e made both for the pop- 
ulation pr operties I Blumenthal et alj[l984l) and for the large-scale 
clustering jDavisetal .11 19851) of galaxies. 

Measurements of the fluctu ation spectru m of the Cosmic Mi- 
crowave Background jSmoot et a l. 1992: Sp ergel et al. 2003 ) and 



of the apparent acceleration of the cosmic expansion jRiess et al.l 
ll998l : |Perlmutter et al.|[l99^ elevated the CDM model, in its vari- 
ant with a cosmological constant (ACDM), to the status of a 
standard paradigm. At the same time improving numerical tech- 
niques and faster computers have enabled detailed simulation of 
the formation and evolution of the galaxy population within this 
paradi gm throughout a sig nificant fraction of the observable Uni- 
verse jspringel et allbOOa) . Nevertheless, direct observational evi- 
dence for halos as extended as the paradigm predicts around galax- 
ies like our own has so far come only from statistical analyses of the 
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dynam ics of satellite galaxies (e.g. lZaritskv et al.ll997l : |prada et al.1 
2003h and of the gravitational lensing of background galaxies (e.g. 
Seliakll2002l : iMandelbaum et al.ll2006l) based on large samples of 
field spirals. 

The earliest observational indication that the effective mass of 
the Milky Way must be much larger than i ts stellar mass came from 
the Timing Argument (hereafter TA) of iKahn & Woltieil ( Il959h . 
These authors noted that the Local Group is dominated by the two 
big spirals, and that these are currently approaching each other at 
about 100 kms^^. (The next most luminous galaxy is M33 which 
is probably about a factor of 10 less massive than M31 and the 
Galaxy.) This reversal of the overall cosmic expansion must have 
been generated by gravitational forces, and since the distance to the 
nearest external bright galaxy is much greater than that between 
M31 and the Milky Way, these forces are presumably dominated 
by material associated with the two spirals themselves. 



iKahn & Woltieil set up a simple model to analyse this situation 
- two point masses on a radial orbit. These were at pericentre (i.e. 
at zero separation) at the Big Bang and must have passed through 
apocentre at least once in order to be approaching today. Clearly 
this requires an apocentric separation larger than the current sepa- 
ration and an orbital period less than twice the current age of the 
Universe. Together these requirements put a lower limit on on the 
total mass of the pair. A more precise estimate of the minimum pos- 
sible mass is obtained from the parametric form of Kepler's laws 
for a zero angular momentum orbit: 
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where r is the current separation, dr/dt is the current relative ve- 
locity, a is the semi-major axis, x is the eccentric anomaly, t is 
the time si nce the Big Bang (the age of the universe) and M is 
total mass jLvnden-B elll 1 1 9 8 If) . Given observationally determined 
values for r, dr/dt and t, these equations have an infinite set of 
discrete solutions for x, i and M labelled by the number of apoc- 
entric passages since the Big Bang. The solution corresponding to 
a single apocentric passage gives the smallest (and only plausible) 
estimate for the mass, which is about 5 x lO^'^ M0 for current es- 
timates of r, dr/dt and t. Note that this is still a lower limit on the 
total mass, even within the simple point-mass binary model, since 
any non-radial motions in the system would increase its present ki- 
netic energy and so increase the mass required to reverse the initial 
expansion and bring the pair to their observ ed separation by the 
present dav (see | Einasto & Lvnden-Belllll983) . 

As lKahn & Woltiet realised, this timing estimate of the total 
mass of the Local group exceeds by more than an order of mag- 
nitude the mass within the visible regions of the galaxies, as es- 
timated from their internal dynamics, in particular, from their ro- 
tation curves. Thus 90% of the mass must lie outside the visible 
galaxies and be associated with little or no detectable light. Mod- 
ern structure formation theories like ACDM predict this mass to be 
in extended dark matter halos with M{r) increasing very roughly 
as r out to the point where the halos of the two galaxies meet. Such 
structures have no well-defined edge, so any definition of their total 
mass is necessarily somewhat arbitrary. In addition, their dynamical 
evolution from the Big Bang until the present is substantially more 
complex than that of a point-mass binary. Thus the mass value re- 
turned by the Timing Argument cannot be interpreted without some 
calibration against consistent dynamical models with extended dark 
halos. 

A first c alibration of this type was carried out by 
iKroeker & Carlberg ( 199 li) using simulations of an Einstein-de Sit- 
ter CDM cosmo gony. Here we use th e very much larger Millen- 
nium Simulation Ispringel et alj2005l) to obtain a more refined cal- 
ibration based on a large ensemble of galaxy pairs with observable 
properties similar to those of the Local Group. We find that the 
standard timing estimate is, in fact, an almost unbiased estimate of 
the sum of the conventionally defined virial masses of the two large 
gala xies. 

IZaritskvetal.l ( ll989[) attempted to measure the halo mass of 
the Milky Way alone by measuring radial velocities for its dwarf 
satellites and assuming the population to be in dynamical equilib- 
rium in the halo potential. They noted, however, that one of the 
most distant dwarfs, Leo L has a very large recession velocity and 
as a result provides a interesting lower limit on the Milky Way's 
mass by a variant of the original Timing Argument. To reach its 
present position and radial velocity, Leo I must have passed peri- 
centre at least once since the Big Bang and now be receding from 
the Galaxy for (at least) the second time. 

Applying the point-mass radial orbit Equations (1) - (3) to 
this case gives a lower bound of about 1.6 x 10^^ Mq. This seems 
likely to be a significant underestimate, since Leo I could not have 
passed through the centre of the Milky Way without being tidally 
destroyed so its orbit cannot be purely radial. Below we calibrate 



the Timing Argument for this case also, finding it to work well 
although with significantly more scatter than for the Local Group 
as a whole. This is because the ACDM paradigm predicts that the 
dynamics on the scale of Leo Fs orbit ('^ 200 kpc) is typically 
more complex than on the scale of the Local Group as a whole 
(~ 700 kpc). 

Our paper is organised as follows. In Section |2] we briefly 
describe the Millennium Simulation and the selection criteria we 
use to define various samples of 'Local Group-like' pairs and of 
'Milky Way-like' halos. In Section[3] we plot the ratio of true total 
mass to Timing Argument mass estimate for these samples, and 
we use its distribution to define an unbiased TA estimator of true 
mass with its associated confidence ranges. In Section [T2I this is 
then applied to the Local Group in order to obtain an estimate its 
true mass with realistic uncertainties. Section [373] attempts to carry 
out a similar calibration for the TA-based estimate of the Galaxy's 
halo mass from the orbit of Leo I. We conclude in Section|4]with a 
summary and brief discussion of our results. 



2 THE MILLENNIUM SIMULATION 

The Millennium Simulation is an extremely large cosmological 
simulation carried out by the Virgo Consortium dSpringel et al.l 
[2005). It followed the motion of iV = 2160^ dark matter par- 
ticles of mass 8.6 x 10* Mq within a cubic box of co- 
moving size 500 Mpc. Its comoving spatial resolution (set 
by the gravitational softening) is 5 h^^ kpc. The simulation 
adopted the concordance ACDM model with parameters Qrn ~ 
0.25, f^fc = 0.045, = 0.73, f^A = 0.75, n = 1 and erg = 
0.9, where, as usual, we define the Hubble constant by Hq = 
lOO/ikms"'^ Mpc~^. The current age of the universe is then 
13.6 x 10^ yr. The positions and velocities of all particles were 
stored at 63 epochs spaced approximately logarithmically in ex- 
pansion factor at early times and at approximately 300 Myr in- 
tervals after z — 2. For each such snapshot a friends-of-friends 
group-finder was used to locate all virialised structures, and their 
self-bound substructu res (subhalos) were identified using SUBFIND 
dSpringel et"aLll200ll) . Halos and subhalos in neighbouring outputs 
were then linked in order to build formation history trees for all the 
subhalos present at each time. These data are publicly available at 
the Millennium release sitfl A "Milky Way" halo at 2 = typi- 
cally contains a few thousand particles and several resolved subha- 
los. 

Galaxy formation was simulated within these merging history 
trees by using semi-analytic models to follow the evolution of the 
baryonic components associated with each halo/subhalo. Processes 
included are radiative cooling of diffuse gas, star formation, the 
growth of supermassive black holes, feedback of energy and heavy 
elements from supemovae and AGN, stellar population evolution, 
galaxy merging and effects due to a reionising UV background. The 
z = galaxy catalogue we analyse here corresponds to the specific 
model of|Croton et al. (200fj) and details of its assumptions and pa- 
rameters can be found in that paper. Data for the galaxy population 
at all redshifts are available at the Millenn ium web site for the up- 
dated model of lDe Lucia & BlaizotI • 2007 ), as well as fo r the inde- 
pendent galaxy formation model of Bower et"ai] ( l2006h . All these 
models are tuned to fit a wide variety of data on the nearby galaxy 
population, and in addition fit many (but not all!) available data at 



http://www.mpa-garching.mpg.de/millennium 



Masses for the Local Group and the Milky Way 3 



higher redshift (see, for example. iKitzbichler & Whit3l2006h . The 
details of the galaxy formation modelling are not, however, impor- 
tant for the dynamical issues which are the focus of our own paper. 

At jz = there are 18.2 x lO'' halos/subhalos identified in the 
simulation to its resolution limit of 20 particles. The galaxy forma- 
tion model populates these with 8, 394, 180 galaxies brighter than 
an absolute magnitude limit of Mb = — 16.7 above which the cata- 
logue can be considered reasonably complete. These catalogues list 
a number of properties for the halos, subhalos and galaxies which 
will be important for us. Galaxies are categorised into three types 
according to the nature of their association with the dark matter. A 
Type galaxy sits at the centre of the dominant or main subhalo and 
can be considered the central galaxy of the halo itself (formally, of 
the FOF group). A Type 1 galaxy sits at the centre of one of the 
smaller non-dominant subhalos associated with a FOF group. Fi- 
nally, a Type 2 galaxy is associated to a specific particle and no 
longer has an associated subhalo because the object within which 
it formed was tidally disrupted after accretion onto a larger halo. 
Such galaxies merge with the central galaxy of their new halo after 
waiting for a dynamical friction time. 

Each galaxy in the catalogue has an associated "rotation veloc- 
ity" Vmax- This is the maximum of the circular velocity Vc(r) = 
{GM{r)/ry^^ of its subhalo for Types and 1; for Type 2 objects 
Vmax is frozen to its value at the latest time when the galaxy still 
occupied a subhalo. Type and 1 galaxies also have an associated 
mass Mhaio which is the mass of the self-bound subhalo which sur- 
rounds them. Finally, halos of Type galaxies have a conventional 
"virial mass" M200, defined as the total mass within the largest 
sphere surrounding them with an enclosed mean density exceeding 
200 times the critical value. Below we will consider both Mhaio 
and M200 as possible definitions for the "true" masses of M3 1 and 
the Galaxy. 

We use the Millennium Simulation to construct samples of 
mock Milky Way/ Andromeda galaxies and of mock Local Groups 
as follows. We begin by identifying all Type or Type 1 galaxies 
with characteristic "rotation velocity" either in the narrow range 
200 ^ Vmax < 250kms~^ or in the wider range, 150 ^ 
Vmax < 300 km s^^. This produces samples of 166,090 and 
699, 177 galaxies respectively. The exclusion of Type 2 galaxies 
reduces the samples by about 5-6% in each case, but the excluded 
galaxies are in any case not plausible analogues for the Local Group 
giants since they are almost all members of large groups or clusters. 
We also consider subsamples in which the morphologies predicted 
by the semianalytic model are forced to approximate those of M3 1 
and the Galaxy. Specifically, we require a bulge-to-total luminos- 
ity ratio in the range 1.2 ^ MB.buige — MB,totai < 2.5 so that 
the disks are 2 to 9 times brighter than the bulges in the B-band. 
This morphology cut reduces the samples in the two Vmax ranges 
to 62, 605 and 271, 857 galaxies respectively. 

We then identify Local Group analogues in each of these 
four samples by identifying isolated pairs with separations in the 
range of 500 — l,000kpc and with negative relative radial ve- 
locities. (Note that this is the true relative velocity rather than the 
relative peculiar velocity, i.e. we have added the Hubble expan- 
sion to the relative peculiar velocity and have required the result 
to be negative.) We identify isolated pairs by keeping only those 
which have no "massive" companion, defined as a galaxy with 
Vmax ^ 150kms~^, within a sphere of 1 Mpc radius centred on 
the mid-point of the binary, and no nearby cluster, defined as a halo 
with M200 > 3 X 10^^ Mq within 3 Mpc of the mid-point of the 
binary. These cuts ensure that the dynamics are dominated by mass 
associated with the two main systems, as appears to be the case 



for the Local Group. For galaxies selected in the narrower Vmax 
range we then find 178 pairs when the morphology cut is applied 
and 1, 128 pairs when it is not. For the wider Vmax range the corre- 
sponding numbers are 2, 815 pairs and 16, 479 pairs respectively. 

When calibrating the TA estimator it proves advantageous to 
use simulated pairs with dynamical state quite close to that of the 
real Local Group. As we will see below, this eliminates some sys- 
tems where the dominant motion is not in the radial direction and 
the TA therefore significantly underestimates the mass. We there- 
fore make one final cut which requires the approach velocity of the 
two galaxies to lie between 0.5 and 1.5 times the value measured 
for the real system (—130 kms^^). This results in our final sets of 
Local Group lookalikes. For the narrower Vmax range we end up 
with 117 pairs when the morphology cut is applied and 758 pairs 
when it is not, while for the wider Vmax range the corresponding 
numbers are 1, 273 pairs and 8, 449 pairs respectively. 

When we study the application of the Timing Argument to 
the Milky Way-Leo I system, we consider individual galaxies from 
both our Vmax ranges. We require these to be isolated by insisting 
that there should be no bright/massive companion (with luminos- 
ity exceeding 10% of that of the host or Vmax > 150kms~^) 
closer than 700 kpc and no massive group (defined as above) closer 
than 3 Mpc. This produces samples of 137, 926 and 266, 229 po- 
tential hosts in the cases with and without the morphology cut for 
the wider Vmax range, and 29, 245 and 57, 816 potential hosts for 
the narrower range. We then search for Leo I analogues around 
these hosts by identifying companions in the separation range 200 
to 300 kpc with Vmax{comp) < 80kms~\ Mb < -16.7 and 
Vra ^ 0.7Vmaa;(host) where Vra is the relative radial velocity 
of the two objects and the last condition reflects the fact that Leo 
I is useful for estimating the Milky Way's mass only because its 
recession velocity is comparable to the Galactic rotation velocity 
(Vra 0.8Vmaii(host) for the real Leo I-Milky Way system). 
Pairs sharing the same MW-like host are excluded in the final list. 

With these cuts we find 344 and 896 satellite-host pairs in 
the samples with and without the morphology cut for the looser 
Vmax range, and 168 and 374 for the tighter range. These relatively 
small numbers reflect the fact that only about 10% of potential hosts 
actually have a faint companion in this distance range which is still 
bright enough to be resolved, and fewer than 5% of these satellites 
are predicted to have positive recession velocities comparable to 
that observed. 



3 RESULTS 

3.1 Calibration of the Timing Argument mass for the Local 
Group 

For each simulated Local Group analogue the separation and rel- 
ative radial velocity of the two galaxies can be combined with the 
age of the Universe (taken to be 13.6 Gyr) to obtain a Timing Argu- 
ment mass estimate Mta (Equations (1) to (3)). The true mass of 
the pair Mt,- is harder to define because of the extended and com- 
plex mass distributions predicted by the ACDM model. The mass 
of an individual dark halo is often taken to be M200 the mass within 
a sphere of mean density 200 times the critical value, so a natural 
choice for Altr is the sum of i\/200 for the two galaxies. The Mil- 
lennium Simulation database only lists M200 for Type galaxies, 
those at the centre of the main subhalo of each friends-of-friends 
particle group. Many of our LG analogues lie within a single FOF 
group. One of the pair is then a Type 1 galaxy, the central object 
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of a subdominant subhalo, and so has no listed value for Khoo- In 
such cases we have gone back to the particle data for the simulation 
in order to measure M200 directly also for these galaxies. 

An alternative convention is to define Mtr as the sum of the 
values of Mhaio, the maximal self-bound mass of each subhalo; 
this is included in the database for both Type and Type 1 galax- 
ies. In the following we use the notation Mtr,200 and Mtr,haio to 
distinguish these two definitions. For either we can calculate the 
ratio of true mass to Timing Argument estimate, 

= Mtr.x/MTA, (4) 

where the suffix x is 200 or halo depending on the definition 
adopted for Mtr- If the Timing Argument is a good estimator of 
true mass, our samples of LG analogues should produce a narrow 
distribution of A values. This distribution then allows the TA mass 
estimate for the real Local Group to be converted into a best esti- 
mate of its true mass, together with associated confidence intervals. 

Our preferred sets of Local Group analogues contain simu- 
lated galaxy pairs which mimic the real system in terms of mor- 
phology, isolation, pair separation and pair approach velocity. In 
addition, they require the halos of the simulated galaxies to have 
Vmax values within about ±10% and ±35% of those estimated for 
M31 and the Galaxy for the tight and loose ranges of Vmax, re- 
spectively. In order to understand the influence of these constraints 
we give results below not only for our "best" samples but also for 
samples where the various constraints are relaxed. Thus, we con- 
sider samples in which 1) both morphology and isolation require- 
ments are applied (our preferred case), 2) the isolation requirement 
is removed, 3) the morphology requirement is removed, and 4) both 
morphology and isolation requirements are removed. For each case, 
we compare results for the two allowed ranges of Vmax and we 
also examine the effect of loosening the radial velocity constraint 

to Vra < 0. 

Fig.[T]gives histograms of the distribution of A200 for a sam- 
ple in the narrow Vmaa; -range with our preferred isolation, mor- 
phology and radial velocity cuts, as well as for three samples with 
the same Vmax and Vra cuts but with reduced morphology and iso- 
lation requirements. Fig.|2]presents these same distributions in cu- 
mulative form and compares them with the corresponding distribu- 
tions for samples with the loosened circular velocity requirement, 
150kms~^ ^ Vmax < 300 km s^^. In both plots black curves 
refer to class (1) samples for which both isolation and morphology 
cuts are imposed, while red, green and blue curves refer to sam- 
ples in classes (2), (3) and (4) respectively. Results for the broader 
Vmax selection are indicated by dashed curves in Fig. [2] We give 
numerical results for various percentile points of these distributions 
in Table [T] and repeat all these in Table |2] for samples where the 
separation velocity requirement has been loosened to Vra < 0. 

The first and most important point to note from these these 
figures and tables is that the median value of A20Q is very robust 
and only varies between 0.98 and 1.34 for our full range of sam- 
ple selection criteria. With our prefened cuts the median values are 
1.15 and 0.99 for the narrow and wide Vmax samples respectively. 
The best estimate of the true mass of the Local Group (for this def- 
inition) is thus very similar to its Timing Argument mass estimate, 
and depends very little on the calibrating sample of simulated pairs. 

The second important point is that the width of the distribution 
of A200 does depend on how the calibrating sample is defined. In 
particular, it is narrower for samples with the more restrictive Vmax 
range, and for given Vmax range it is smallest for samples with our 
preferred cuts, those which match the dynamical and morpholog- 
ical properties of the Local Group most closely. For the narrow 
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Figure 1. Normalised histograms of A200, the ratio of true mass to 
Timing Argument estimate, for samples of Local Group analogues with 
200kms-i ^ Vmax < 250kms-i and -195kms-i < Vra < 
—65 kms^^. The black histogram refers to our preferred selection where 
both isolation and morphology requirements are imposed. For the red his- 
togram the isolation requirement has been removed, for the green histogram 
the morphology requirement, and for the blue histogram both requirements. 
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Figure 2. Cumulative distributions of j4200 the ratio of true mass to Tim- 
ing Argument estimate for Local Group analogues with —195 km s^^ < 
Vra < — 65kms^^. The solid curves correspond to the four samples 
already plotted in Fig. [T] while the dashed curves are for samples with 
150kms~t ^ Vmax < 300 km s"^. The colour coding is the same as in 
Fig.[T] black indicates samples with our preferred isolation and moiphology 
constraints. 
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Vmax sample the interquartile range is a factor of just 1.6, and the 
upper and lower 5% points are separated by a factor of 3. For the 
wider velocity range the interquartile range is a factor of 1.8 and 
the 5% points are separated by a factor of 5.7. This shows the Tim- 
ing Argument to be remarkably precise for systems similar to the 
Local Group. 

The broadening of the distribution as the selection require- 
ments are relaxed is easy to understand. Removing the isolation 
requirements allows third bodies to play a significant role in the 
dynamics. This can extend the upper tail of the ^200 distribution if 
mass from the third body falls inside -R200 for one of the pair galax- 
ies or if the gravity of the third galaxy produces a tidal field which 
opposes the attraction between the pair members. It can extend the 
lower tail if the mass of the third body lies between the pair mem- 
bers but outside their i?200 spheres, thus enhancing their mutual 
attraction without adding to their mass. Removing the morphology 
constraint moves the whole distribution towards larger values and 
this effect is most pronounced in the large A200 tail. This is because 
objects with more dominant bulges have more complex merger his- 
tories. They typically form in denser regions and their halos tend to 
be more massive and to have more complex structure. 

Loosening the requirements on Vmax affects the distribution 
in a complex way. There is a tight relation between Vmax and 
M200 (also Mhaio) in the AC DM structure formation model (e.g. 
iNavarro. Frenk &" Whitelll997h . Thus if we place tight restrictions 
on the Vmax values of our galaxies, we will obtain a sample of Lo- 
cal Group analogues with a nan^ow range of Mtr values. If, in addi- 
tion, we force the parameters which enter in the Timing Argument 
(the pair separation and relative radial velocity) to lie in narrow 
ranges, then the TA mass estimate itself is tightly constrained. The 
distribution of yl2oo is thus forced to be narrow as a consequence 
of our selection criteria. 

A second effect is that most of the new pairs added by 
widening the requirement on Vmax have at least one galaxy with 
150 kms^^ ^ Vmax < 200 kms^^, thus with relatively low A/tr. 
This simply reflects the strong dependence of halo abundance on 
Vmax- Given that halo mass scales approximately as Vma,x it is 
striking that the addition of so many pairs containing a "low mass" 
galaxy reduces the median value of A200 by just 15%. The low tail 
of the distribution is more strongly affected, by almost a factor of 
2 at the lower 5% point, but the upper end of the distribution is 
barely affected at all. This demonstrates that the main body of the 
distribution is weakly affected by restrictions on Vmax, but that the 
lower tail (which is needed to place a lower limit on the true mass of 
the Milky Way) is suppressed if Vmax is not allowed to take small 
values. 

In Table |2] we show the effect of weakening the cut on rela- 
tive radial velocity to require only that the two main galaxies be 
approaching. Again this has remarkably little effect on the median 
A200 values. A comparison with Table [T] shows them all to be in- 
creased by about 10%- 15%. The effects on the tails of the distribu- 
tions are more substantial. The 95% point is typically increased by 
about a factor of 2. This is because the sample now includes a sub- 
stantial number of pairs with small Vra (and thus smaller TA mass 
estimate) for which tangential motion is important for their current 
orbit. The 5% point of the distribution is significantly reduced, re- 
flecting the fact that our restrictions on relative approach velocity 
exclude a non negligible number of systems with approach veloci- 
ties larger than 195 km s~^, and thus with large TA mass estimates. 
Such systems must have more mass outside the conventional virial 
radii of the two galaxies than do typical Local Group analogues in 
our samples. 



In conclusion, we believe our most precise and robust estimate 
of the distribution of A200 to be that obtained with our preferred 
morphology, isolation and radial velocity cuts for 150 km s^^ ^ 
Vmax < 300 kms^^, and we will use this distribution in the next 
section to estimate the true mass of the Local Group. Although the 
tails of the distribution are suppressed still further for a narrower 
range of Vmax, this is at least in part due to the artificial effects 
mentioned above. In addition the number of Local Group analogues 
is too small in this case for the tails of the distribution to be reliably 
determined. From Table [T] we see that the best estimate of the true 
mass of the Local Group (which we take to be that obtained using 
the median value of A200) is almost identical to the direct TA esti- 
mate. The most probable range of true mass (given by the quartiles 
of A200) extends to values about 30% above and below this, while 
the 95% confidence lower limit on the true mass (given by the 5% 
point of the ^200 distribution) is a factor of 2.9 smaller. 



3.2 Application to the Local Group 

The three observational parameters needed to make a Timing Ar- 
gument mass estimate for the Local Group are the separation be- 
tween the two main galaxies, their radial velocity of approach and 
the age of the Universe. The latter is now determined to high preci- 
sion through measur ements of microwave background fluctuations. 
ISpergel et all i IiOOTI) give 13.73 ± 0.16 Gyr. The distance to M31 
is also kn own to high precision . We a dopt the value 784 ±21 kpc 
given bv IStanek & GarnavichI h998l) based on red clump stars, 
noting that it agrees al most exactly w ith the slightly less pre- 
cise value obtained by iHollanj l ll998l) from fits to the colour- 
magnitude diagrams of M31 globular clusters. Although the he- 
liocentric recession velocity of M3 1 is known even more precisely 
(—301 ± 1 kms^^ according to lCourteau & van den Bergrill99a) 
the approach velocity of the two giant galaxies is less certain be- 
cause of the relatively p oorly known rotation velocity of the M ilky 
Way at the Solar radius. Ivan der Marel & Guhathakurtal | |2007|) go 
through a careful analysis of the uncertainties and conclude that 
Vra = 130 ± 8kms~^. Inserting these modem values into Equa- 
tions (1) to (3) we obtain our Timing Argument estimate of the 
mass of the Local Group: 

Mlg,ta ~ 5.32 ± 0.48 x 10^^ M© , (5) 

where the uncertainty is dominated by that in the relative radial 
velocity. This uncertainty is still small in comparison to the scatter 
in the ratio of true mass to TA estimate, so we will neglect it in the 
following. The apocentric distance of the implied relative orbit of 
the two galaxies is 1103 ± 30 kpc. 

We now combine this Timing Argument estimate with the dis- 
tribution of ^200 obtained in the last section for our most precise 
and reliable sample of Local Group analogues (the sample with our 
preferred morphology, isolation and radial velocity cuts and with 
the wider allowed range of Vmax) to obtain our best estimate of the 
true mass of the Local Group, defined here as the sum of the M200 
values of the two main galaxies: 

MLCtrue = 5.27 X 10'^ Mo , (6) 

or log MLG,true/MQ — 12.72. The most plausible range for this 
quantity is then [12.58, 12.83] with a 95% confidence lower limit 
of 12.26, i.e MLG,true > 1.81 x 10^^ Mq at 95% confidence. 
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Table 1. Percentage points of the A200 distribution for samples of LG analogues with 
-195kms-i < Vra < -65kms-i 





5% 


25% 


50% 


75% 


95% 


# of pairs 


200kms-l ^ V^ax < 250kms-i 














morphology, isolation 


0.67 


0.93 


1.15 


1.47 


2.05 


117 


morphology, no isolation 


0.61 


0.93 


1.14 


1.52 


2.09 


155 


no morphology, isolation 


0.67 


0.97 


1.20 


1.50 


2.32 


758 


no morphology, no isolation 


0.63 


0.96 


1.22 


1.55 


2.54 


1015 


150kms~l < Vmax < SOOkms^i 














morphology, isolation 


0.34 


0.72 


0.99 


1.27 


1.93 


1273 


morphology, no isolation 


0.33 


0.68 


0.98 


1.29 


2.00 


1650 


no morphology, isolation 


0.41 


0.81 


1.09 


1.40 


2.21 


8449 


no morphology, no isolation 


0.35 


0.77 


1.08 


1.43 


2.41 


11838 



Table 2. Percentage points of the A200 distribution for samples of LG analogues with Vra < 





5% 


25% 


50% 


75% 


95% 


# of pairs 


200kms-l < Vmax < 250kms-i 














morphology, isolation 


0.54 


0.97 


1.33 


1.66 


3.93 


178 


morphology, no isolation 


0.45 


0.94 


1.26 


1.66 


3.93 


241 


no morphology, isolation 


0.54 


1.01 


1.34 


1.82 


4.62 


1128 


no morphology, no isolation 


0.42 


0.96 


1.34 


1.93 


5.11 


1596 


150kms-l < Vmax < 300kms-i 














morphology, isolation 


0.28 


0.85 


1.19 


1.64 


3.30 


2815 


morphology, no isolation 


0.22 


0.77 


1.16 


1.64 


3.38 


3532 


no morphology, isolation 


0.31 


0.89 


1.23 


1.76 


4.06 


16479 


no morphology, no isolation 


0.18 


0.79 


1.19 


1.78 


4.46 


23429 



3.3 Application to the Milky Way 

We now calibrate the IZaritskv et al. I ( Il989h Timing Argument 
which estimates the mass of the Milky Way from the position and 
velocity of Leo I. This again assumes a radial Keplerian orbit, but 
Leo I is taken to have passed through pericentre and to be currently 
moving towards apogalacticon for the second time. Equations (1) to 
(3) then give a unique mass estimate Mmw,ta for the system for 
any assumed distance and radial velocity. This is taken as the Milky 
Way's mass since the mass of Leo I is neghgible in comparison. 

Proceeding as for the Local Group, we select Milky Way - 
Leo I analogues from the Millennium Simulation in order to study 
the relation of this TA estimate to the true mass of the Milky Way, 
which we again take to be Af200- Thus we define the ratio 

B200 ~ M200 / Mmw,ta (7) 

and investigate its distribution in various samples of analogue host- 
satellite systems. In particular, we consider samples of isolated host 
galaxies (as defined in Section[2]l using both our looser and tighter 
Vmax ranges, both with and without cuts on central galaxy mor- 
phology, and requiring the distance, radial velocity and maximum 
circular velocity of the satellite to satisfy 200 kpc < r < 300 kpc, 
Vra > 0.714iaa:(host) and Vmax{comp) < 80kms"^. 

Results for these four samples are given in Table|3]and the cor- 
responding cumulative distributions of B200 are plotted in Fig. [3] 
Scatter plots of M200 against Mmw,ta for the four samples are 
shown in Fig.|4l The behaviour is quite similar to that of the Local 
Group TA mass estimator yl200- The median value of B200 is ro- 
bust and varies very little as the definition of the analogue sample is 
changed. Again it is 10 - 15% smaller for samples with the looser 
Vmax selection. Unlike the Local Group case, the median value 



of B20Q is about 1.6 and so is significantly larger than unity. This 
shows that Mmw,ta is biased low as an estimator of true Milky 
Way mass, reflecting the fact that tangential motions are often sig- 
nificant for satellites at the distance of Leo L Assuming a purely 
radial orbit then results in an underestimate of the mass. 

The width of the distribution of B200 is significantly greater 
for samples with the looser Vmax selection, primarily through an 
extension of the tail towards low values. This resembles the be- 
haviour we saw above for A200 but it must have a different cause, 
since our selection criteria for Milky Way analogues put no upper 
bound on Vra, instead placing a lower limit on Vra /Vmax- As a 
result they do not force an upper limit on Mmw,ta of the kind 
imposed on Mlg,ta by our upper limit on Vra for Local Group 
analogues. Fig. |4] shows that the tail of low i3200 values for the 
wider Vmax range is caused a relatively small number of systems 
for which Mmw,ta is anomalously large. These are objects with 
anomalously large values of Vra and seem to occur preferentially at 
small M200, corresponding to values of Vmax below 200 kms^^. 

The bulk of the points in Fig. |4] scatter fairly symmetrically 
about the median relation M200 = l-6Afi\/vK,TA which we show 
as a dashed straight line. Their mean slope is somewhat steeper than 
strict proportionality because our distance constraint on "Leo Fs" is 
expressed in units of kpc rather than of R200 or Vmax/ Hq. Distant 
outliers occur only the low M200 side of this relation, suggesting 
that they may be a consequence of resolution problems in the Mil- 
lennium Simulation. For Vmax ~ 150 km s^^, typical halos are 
represented by fewer than 1, 000 particles and it seems likely that 
difficulties in describing the dynamics of their satellite substruc- 
tures may begin to surface. In addition, the sample sizes are rela- 
tively small, particularly when we impose a morphology cut, so that 
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Table 3. Percentage points of the -B200 distribution for samples of MW - Leo I analogues with 

Vra ^ O.TVmaiChost) 





5% 


25% 


50% 


75%) 


95% 


# of pairs 


200kms-l < V^ax(^ost) < 250 km s"! 














morphology 


0.71 


1.27 


1.71 


2.01 


2.62 


168 


no morphology 


0.71 


1.25 


1.67 


2.01 


2.55 


374 


150kms-l ^ Vmaxihost) < 300 km s"! 














morphology 


0.39 


1.04 


1.50 


1.89 


2.47 


344 


no morphology 


0.51 


1.14 


1.55 


1.98 


2.66 


896 



the estimates of the tails of the distributions may be noisy. This may 
explain in part the apparent excess of outliers in the morphology- 
selected sample with the wider Vmax range. 

The observational data needed to obtain the TA estimate of 
the Milky Way's mass are the age of the Universe and the Galac- 
tocentric distance and radial velocity of Leo I. As a bove, we take 
the ag e of the Universe to be 13.73 ± 0.16 Gyr from lSpergel et al.] 
( I2OO7 ). For the heliocentric d istance to Leo I we adopt 254± 19 kpc 
from Bellazzini et alj | |2004|) . The heliocentric radial velocity of 
Leo I is very precisely determined, 283 ±0.5 kms^^ according to 
iMateo. Olszewski & Walked ( I2OO7I) . Based on an assumed Galactic 
rotation speed at the Sun of 220 ± 15kms^^, we derive a corre- 
sponding Galactocentric radial velocity of 175 ± 8kms^^. When 
substituted into Equations (1) to (3), these parameters produce a TA 
estimate for the Milky Way's mass of 



LOO 



0.75 



0.50 



0.25 




Mmw,ta = 1.57 ± 0.20 X lO'^ Mq . (8) 

As was the case for the Local Group, the fundamental observational 
quantities are so well defined that the uncertainty of this estimate is 
much smaller than the expected scatter in B200 • We will therefore 
neglect it in the following. The implied apocentric distance of Leo 
I is 619 ± 26 kpc. Since this is about half the apocentric distance of 
the M3 1 - Milky Way relative orbit in the TA model of SectionlXIl 
perturbations of the orbit of Leo I due to the larger scale dynamics 
of the Local Group seem quite likely. 

For the reasons discussed above, we consider our most precise 
and robust estimate for the distribution of ^200 to be that obtained 
for host galaxies with 150 kms^^ ^ V,nax{iiost) < 300 kms^^ 
and with no morphology cut. The median of this distribution then 
gives our best estimate of the true halo mass of the Milky Way: 



nonl I -T I 

0.1 1.0 10.0 

1^200 

Figure 3. Cumulative distributions of B200 the ratio of true Milky Way 
mass (taken to be Af20o) to TA estimate for four samples of isolated Milky 
Way - Leo 1 analogues from the Millennium Simulation. The red curve 
refers to Milky Way analogues with 200kms~^ ^ Vmax {host) < 
250 km and with morphology matching the Milky Way. For the black 
curve the circular velocity requirement is loosened to 150kms~^ ^ 
Vmax (host) < 300kms~^, for the blue curve the morphology require- 
ment is removed, and for the green curve both requirements are relaxed. In 
all cases we require Vra ^ 0.7Vmaa; (host). 



M2oo,MW = 2.43 X 10'^ Mq , (9) 

or log M200 /Mq = 12.39. The quartiles of the distribution imply 
[12.25, 12.49] for the most probable range of this quantity, while 
the 5% point implies a lower limit of 11.90 at 95% confidence. 
Thus the implied mass of the Milky Way is roughly half that of the 
Local Group as a whole, as might be expected on the basis of the 
similarity of the two giant galaxies. It is quite similar to other recent 
estimates based on applying equilibrium dyna mics to the system of 
distant Milky Way satellites and halo st ars (e.g. lWilkinson & Evanj 
ll999l : ISakamoto. Chiba & Beersll2003h . A significantly smaller es- 
timate came from the a nalysis of the high -velocity tail of the lo- 
cal stellar population bv lSmithetai]j2007h . but we note that such 
analyses, in reality, only place a lower limit on the mass of the 
halo, since the distribution of solar neighbourhood stars may well 
be truncated at energies significantly below the escape energy. 



3.4 An alternative mass measure? 

The halo masses we have quoted so far have been based on the 
"virial masses" M200 of simulated halos. This choice is, of course, 
somewhat arbitrary, and it may not correspond particularly well to 
the radii within which individual isolated halos are approximately 
in static equilibrium. As an alternative convention, we here con- 
sider defining the mass of an individual halo to be that of the corre- 
sp onding self-bound sub halo identified by the SUBFIND algorithm 
of ISpringel et al.l ( l200lh . This algorithm typically includes mate- 
rial outside the radius R200 within which M200 is measured, but it 
excludes any material which is identified as part of a smaller sub- 
halo orbiting within the larger system. In this paper we denote this 
subhalo mass as M^aio- 

In the left panel of Fig. |5] we plot Mhaio against Af200 for all 
halos in our preferred sample of Local Group analogues, that with 
our preferred morphology, isolation and radial velocity cuts and 
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Figure 4. Scatter plots of the TA estimate of the "Milky Way's" mass against its true mass, M2oo- The left panel is for isolated host galaxies with morphology 
matching that of the real Milky Way, while the right panel is for samples with no such morphology restriction. In both panels red dots refer to simulated 
host galaxies with 200 km ^ Vma^Chost) < 250 km s~^, while black dots indicate other hosts in the broader range 150 km ^VmaxQ^ost) < 
300 kms"^ In all cases the analogue of Leo I is required to have 200 kpc < r < 300 kpc, Vl-a ^ 0.7Knaa;(host) and Vmaa;{comp) ^SOkms ^.The 
straight dashed lines in the two panels show the approximate median relation M200 = l-6Mjv./vK TA 



with 150 kms~^ ^ Vmax < 300 km s~^. Black and red points in 
this plot refer to Type and Type 1 subhalos respectively. The right 
panel of Fig. |5]is a similar plot for the "Milky Way" halos in our 
preferred sample of Milky Way - Leo I analogues, again the sam- 
ple which is matched in morphology and which has the wider Vmax 
range. In both panels it is clear that the correspondence between the 
two mass definitions is quite tight, and that Mhaio tends to be some- 
what larger than i\/2oo • In addition the left panel shows that Type 1 
halos have smaller Mhaio for given A/200 than do Type halos, as 
would naively be expected. The average value of log Mhaio/M2oo 
for the halos in the left panel is 0.100 for the Type O's and —0.004 
for the Type I's, while it is 0.079 for the "Milky Way" halos in the 
right panel. 

This close correspondence between the two mass definitions 
carries over to the distribution of our ratios of "true" mass to timing 
mass. In Table |4] we give percentage points of the Ahaio distribu- 
tion for the 8 samples of Local Group analogues already consid- 
ered above. They can be compared directly with the numbers given 
in Table [T] for these samples. To a good approximation the distri- 
bution of Ahaio agrees with that of j42oo except that all values are 
shifted upwards by about 16-20%. 

The same is also true for estimates of the Milky Way's mass 
obtained using the TA applied to Leo I. This can be seen from Ta- 
ble [5] which repeats Table [3] except that we now give percentage 
points for Bhaio rather than £200 • Clearly it is of rather little impor- 
tance which definition of halo mass we adopt: the results obtained 
with our two definitions are very similar. 



4 DISCUSSION AND CONCLUSIONS 

The statistical argument underlying the analysis of this paper is 
more subtle than it may at first appear, so it is worth restating it 



somewhat more formally in order to understand what is being as- 
sumed in deriving the mass estimates for the Local Group and for 
the Milky Way given above. 

We believe that the mass distributions around galaxies are 
much more extended than the visible stellar distributions, and that 
these have been assembled from near-uniform "initial" conditions 
in a manner at least qualitatively resembling that in a ACDM 
universe. Thus the assembly histories of the Local Group and of 
the Milky Way's halo differ in majo r ways from tho s e assu med 
by the original Timi ng Arguments of iKahn & Woltieil jl959l) and 
IZaritskveta"i] ( ll989h . In addition, the meaning of the derived mass 
values needs clarification. We wish to use the Millennium Simula- 
tion to calibrate the TA estimates against conventional measures of 
halo mass, and to test the general applicability of the Timing Ar- 
gument. However, we want to do this in a way which avoids any 
significant dependence on the details of the ACDM model, for ex- 
ample, on the exact density profiles, abundances and substructure 
properties which it predicts for halos. 

Our method uses the simulation to estimate the distribution 
of the ratio of "true" mass to TA mass estimate for samples of 
objects whose properties "resemble" those of the observed Local 
Group and Milky Way - Leo I systems. Our restrictions on sepa- 
ration and radial velocity implement this similarity requirement in 
a straightforward way, but our constraints on Vmax have a more 
complex effect. Although the true Vmax values for M31 and the 
Milky Way are very likely within our looser range (150 kms"^ ^ 
Vmax < 300kms~^) the simulation exhibits a tight correlation 
between Vmax and M200 ■ Imposing fixed limits on Vmax is thus 
effectively equivalent to choosing a fixed range of M200 • As a re- 
sult, we are in practice estimating the distribution of A200 or -6200 
for systems of given true mass, subject to the assumed constraints 
on separation and radial velocity. However, when we apply our re- 
sults to estimate true masses for the Local Group and the Milky 
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Figure 5. In the left panel we plot Mh^lo against M200 for all the galaxy (sub)halos in our sample of Local Group analogues with our preferred moiphology, 
isolation and radial velocity cuts and with 150kms~^ ^ Vmax < 300kms~^. Type subhalos are plotted black while Type 1 subhalos are red. The 
diagonal line is the one-to-one relation. The right panel is a similar plot in for the "Milky Way" halos in our preferred sample of Milky Way - Leo I analogues. 



Table 4. Percentage points of the Aj-^^i^ distribution for samples of LG analogues with 
-195kms-i < Vra < -65kms-l 





5% 


25% 


50% 


75% 


95% 


# of pairs 


200kms-i < Vniax < 250kms-i 














morphology, isolation 


0.78 


1.11 


1.36 


1.72 


2.38 


117 


morphology, no isolation 


0.69 


1.14 


1.36 


1.77 


2.39 


155 


no morphology, isolation 


0.76 


1.17 


1.44 


1.75 


2.74 


758 


no morphology, no isolation 


0.71 


1.15 


1.44 


1.82 


2.97 


1015 


150kms-l < Vmax < 300kms-i 














morphology, isolation 


0.41 


0.88 


1.20 


1.53 


2.26 


1273 


morphology, no isolation 


0.38 


0.84 


1.19 


1.54 


2.38 


1650 


no morphology, isolation 


0.49 


0.99 


1.31 


1.67 


2.60 


8449 


no morphology, no isolation 


0.37 


0.92 


1.30 


1.71 


2.85 


11838 



Way, we implicitly assume that our distributions of ^200 and B200 
are appropriate for samples of given TA mass estimate, again sub- 
ject to our constraints on separation and radial velocity. It is thus 
important to understand when these two distributions can be con- 
sidered the same. 

The relation can be clarified as follows. From the simula- 
tion we compile the distribution of Mtr /Mta, or equivalently of 
A = In Altr — In Mta, for systems with In Altr in a given range. 
We then implicitly assume that this distribution does not depend on 
Mtr, at least over this range, so that the result can be taken as an 
estimate of the probability density of A at given Altr - Bayes The- 
orem then gives us the probability density function (pdf) for A at 
fixed A'Ita '■ 



/[A I In Mta] 



/[A, In Mta] 
/[In Mta] 
/[A,lnMt,] 
/[In Mta] 



/[A|lnMt.]/[lnMt.] 
f [In Mta] 
= /[A|lnA/t,] 



(10) 



The first line here simply writes the conditional pdf of A at given 
AIta in terms of the joint pdf of the two quantities and the pdf of 
A/ta- The second line then rewrites the joint pdf in terms of the 
equivalent variables A and Mtr, using the fact that the Jacobian of 
the transformation is unity. The third line re-expresses the joint pdf 
as the product of the pdf of A at given Altr times the pdf of Mtr- 
The final line then follows from the normalisation condition, pro- 
vided that / [In Mtr] is constant and / [A | In Altr] is independent of 
Altr - Thus, when estimating Mtr from Mta, we assume a uniform 
prior on In Altr and that the distribution of A does not depend on 
true mass. Both these assumptions appear natural and appropriate. 

The analysis underlying the Timing Argument (Equations (1) 
- (3)) assumes that the relative orbit of the two objec ts is bound 
and ha s conserved energy since the Big Bang. Recentlv. lSales et al] 
( I2OO7I) have shown that in ACDM models this assumption is signif- 
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Table 5. Percentage points of the Bf^^i^ distribution for samples of MW - Leo 1 analogues with 

Vra ^ O.TVmax.MW 





5% 


25% 


50% 


75% 


95% 


# of pairs 


200kms-l < V^nax < 250kms-i 














morphology 


0.90 


1.53 


1.98 


2.40 


3.14 


168 


no morphology 


0.90 


1.50 


1.97 


2.39 


3.21 


374 


150kms-i < Vmax < 300kms-i 














morphology 


0.42 


1.25 


1.81 


2.31 


3.09 


344 


no morphology 


0.61 


1.37 


1.88 


2.36 


3.23 


896 



icantly violated for a substantial number of satellites within halos 
comparable to that of the Milky Way. In particular, they demon- 
strate the presence of a tail of unbound objects which are being 
ejected from halos as a result of 3-body "slingshot" effects during 
their first pericentric passage. These objects are typi cally receding 
rapidly from their "Milky Way", as assumed in the IZaritskv et al.l 
( Il989h argument, but they violate its assumption that the present 
orbital energy can be used to infer the period of the initial orbit (i.e. 
the time from the Big Bang to first pericentric passage). Clearly 
such objects should also be present in the Millennium Simulation, 
although lack of resolution might make t hem under-represen ted in 
comparison to the simulations analysed bv lSales et alj ( l2007l) . Thus 
our analysis takes the possibility of such ejected satellites into ac- 
count, at least in principle. Objects of this type will show up as sys- 
tems with anomalously large TA estimates for their halo mass, and 
Fig.|4]shows a number of outliers which could well be explained in 
this way. Issues of this kind do not effect TA-based estimates of the 
mass of the Local Group since the two big galaxies are currently 
approaching for the first time. 

The only kinematic information about the relative orbit of 
M31 and the Mi lky Way used in ou r analysis is their current ap- 
proach velocity. Ivan der Marel & Gu hathakurtj ( |2007|) show that 
geometric arguments can already constrain the transverse compo- 
nent also, and future astrometry missions such as SIM might allow 
Vtr to be measured directly. Thus it is interesting to ask if our TA 
mass estimate could be significantly refined by measuring the full 
3-D relative motion of the two galaxies, rather than just its radial 
component. We address this in Fig. [6] which plots A200, the ra- 
tio of true mass to TA estimate, against Vtr for a sample of Local 
Group analogues with our preferred morphology, isolation and ra- 
dial velocity cuts, and with 150kms~^ Vmax < 300 kms^^. 
The median Vtr for this samp l e is 8 6kms~^, comparable to the 
Ivan der Marel & Guhathakurt^ ^20071) estimate for the real system. 
There is no apparent correlation of A200 with Vtr, and indeed, the 
medians of A200 for the high and low Vtr halves of the sample are 
both close to 1 and do not differ significantly. Pairs with high Vtr 
do show larger scatter in A200 than pairs on near-radial orbits. For 
given separation, radial velocity and age, the Kepler model implies 
a mass which increases monotonically with Vtr- The absence of 
a detectable trend in Fig. [6] shows that uncertainties in Vtr do not 
dominate the scatter in our TA mass estimate, and that a measure- 
ment of Vtr will not substantially increase the precision with which 
the true mass can be measured. 

In conclusion, our analysis shows the Timing Argument to 
produce robust estimates of true mass both for the Local Group and 
for the Milky Way, as long as "true mass" is understood to mean 
the sum of the conventional masses of the major halos. For the 
Local Group as a whole, the estimate and confidence limits given 
in Section [T2I and in the Abstract appear reliable given the excel- 




100 150 200 
IV„ I [km s-'] 
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Figure 6. Scatter plot of A200 versus transverse velocity for Local Group 
analogues in our sample with prefen'ed morphology, isolation and radial 
velocity cuts and with 150 km ^ 5C Vmax < 300 km s~ ^ . The vertical 
dashed line indicates the median Vtr - The distributions on either side of 
this line are each further split in half at the median values of Vtr (the solid 
horizontal lines.) There is essentially no correlation in this plot, indicating 
that a measurement of the transverse velocity will not significantly improve 
the TA mass estimate. 



lent statistics provided by the Millennium Simulation, the lack of 
any substantial dependence on our isolation and morphology cuts, 
and the relatively simple dynamical situation. Although the results 
based on Leo Fs orbit also appear statistically sound, the more 
complex dynamical situation offers greater scope for uncertainty, 
particularly when trying to place a lower limit on the mass of the 
Milky Way's halo. On the other hand, our best estimate of this mass 
is just under half of our estimate of the sum of the halo masses of 
M3 1 and the Galaxy. Thus the picture presented by the data appears 
quite consistent, and gives no reason to be suspicious of the Milky 
Way results. 
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